Accepted for publication in ApJL. 

Preprint typeset using 1^1^^ style emulateapj v. 08/13/06 



O 

(N 

< 

(N 

o 

u 

6 



(N 
> 

oo 

p 

rn 
O 

o 



ONE MOMENT IN TIME - MODELING STAR FORMATION IN THE ANTENNAE 

Simon J. Karl\ Thorsten Naab^'^, Peter H. Johansson^, Hanna Kotarba^, Christian M. Boily^, Florent 

Renaud^'", Christian Theis"''^ 

^ Universitats-Sternwarte Miinchen, Scheinerstr. 1, D-81679 Miinchen, Germany; 

skarlOusm . Imu . de 

^ Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild-Str. 1, D-85741 Garching bei Miinchen, Germany; 

naabOmpa-gar ching . mpg . de 
^ Observatoire astronomiquc, Univcrsite de Strasbourg and CNRS UMR 7550, 
11 rue de I'Universite, F-67000 Strasbourg, France 
*Institut fiir Astronomic der Univ. Wien, Tiirkenschanzstr. 17, A-1180 Vienna, Austria and 
^ Planetarium Mannheim, Wilhelm-Varnholt- Alice 1, D-68165 Mannheim, Germany 
Accepted for publication in ApJL. 

ABSTRACT 

We present a new high-resolution N-body/SPH simulation of an encounter of two gas-rich disk 
galaxies which closely matches the morphology and kinematics of the interacting Antennae galaxies 
(NGC 4038/39). The simulation includes radiative cooling, star formation and feedback from SNII. 
The large-scale morphology and kinematics are determined by the internal structure and the orbit of 
the progenitor disks. The properties of the central region, in particular the starburst in the overlap 
region, only match the observations for a very short time interval (Ai « 20 Myr) after the second 
encounter. This indicates that the Antennae galaxies are in a special phase only about 40 Myr after 
the second encounter and 50 Myr before their final collision. This is the only phase in the simulation 
when a gas-rich overlap region between the nuclei is forming accompanied by enhanced star formation. 
The star formation rate as well as the recent star formation history in the central region agree well with 
observational estimates. For the first time this new model explains the distributed extra-nuclear star 
formation in the Antennae galaxies as a consequence of the recent second encounter. The proposed 
model predicts that the Antennae are in a later merger stage than the Mice (NGC 4676) and would 
therefore lose their first place in the classical Toomre sequence. 

Subject headings: galaxies: evolution — galaxies: individual (NGC 4038/39) — galaxies: interactions 
— galaxies: star formation — methods: numerical 



1. INTRODUCTION 

In the local universe (z < 0.3) about ~ 5 — 10% of 
all galaxies are interacti ng and merging (e.g. iLotz et al] 
120081 : iBridge et al.ll2010[ ). Mass assembly via this mech- 
anism was more important at earlier cosmi c times when 
majo r mergers were more f requent (e.g. iPatton et al.l 
120021 : IConseh ce et^ I2003D and also more gas-rich 
(e.g. iTacconi e t al.l[2010l ). Major mergers dramatically 
affect the formation and evolution of galaxies. By 
inducing tidal torques they can effi ciently transport gas 
to th e centers of the g alaxies (|Barnes fc Hernquls^ 
119961: iNaa b et all |2006|) . trigger star formation 
(iMihos fc^ crnguis^ 
2008') , fee d super-ir i assive 



iSoringell [2000t iCox et al 
black holes jHopk ins et al 



2005; S pringel et al.l 120051 : [Johansson el a l. 20091) and 
convert spira l galaxies into intermediate-mass ellipticals 
(iBarneijIi gga INaab fc Burkertj|2003l:lRothberg fc Josephl 
[2004; Naa b fc Ostrikeill2009| ) . 

The Antennae galaxies (NGC 4038/39) are the nearest 
and best-studied example of an on-going major merger 
of two gas-rich spiral galaxies. The system sports a 
beautiful pair of elongated tidal tails extending to a pro- 
jected size of ~ 20' (i.e. 106 kpc at an assumed distance 
of 22 Mpc), together with two clearly visible, still dis- 
tinct galactic disks. The latter has been assumed to be 
an indication of an early me rger stat e , putt ing the sys- 
tem in the first place of the iToomrd ()1977[ ) merger se- 
quence of 11 prototypical mergers. Due to their prox- 
imity and the ample number of high-quality observa- 



tions covering the spectrum from radio to X-ray (e.g. 
, Neff fc Ulvcs tad 2000: Wa ng et al..2004: Whitmore et all 
ll999l:lHibbard et al.l 120051: IZezas et 311120061 ) the Anten- 
nae provide an ideal laboratory for understanding the 
physics of merger-induced starbursts through compari- 
son with high-resolution simulations. 

At the center of the Antennae galaxies, HST imag- 
ing has revealed a large number of bright young star 
clusters (> 1000) which plausibly have formed in sev- 
eral bursts of star form ation induced by the interaction 
(jWhitmore et al.l I1999D . The spatial distribution and 
the age of these clusters are correlated: the youngest 
clusters are found in the overlap region (r < 5 Myr), 
while the young starburst is generally located in the over- 
lap and a ring-like configuration in the disk of NGC 
4038 (r < 30 Myr). An intermediate-age population 
(r — 500 — 60 Myr) is distributed t hroughout the dis k 
of NGC 4038 (iWhitmore et al.lll99llZhang et al.ll200l . 

Of particular interest is the spectacular nature of an 
extra-nuclear starburst observed in the du sty overlap re- 
gion between the merging galactic disks (jMirabel et al.l 
|1998; Wang et al. 2004). The Antennae seem to be the 
only interacting system where an off-center st arburst is 
outsh ining the galactic nuclei in the mid-IR ()Xu et al.l 
l2000r ) and among only a few systems w hich show en- 
hanc ed inter-nuclear gas concentrations (jTacconi et al.l 
Il999f ). To date, this prominent feature has not been re- 
produced in any simulat ion of the Antennae system (see 
iBarnes fc HibbardI |2Q09| ). Thus, the question remains 
whether this feature cannot be captured by current sub- 
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Fig. 1. — Simulated gas properties projected on top of Hi kinematic data bv lHibbard et aTl 1 1200 II) at the time of best match. Simulated 
gas particles are displayed in blue (NGC 4038) and red (NGC 4039). Yellow points represent the observational data. Upper left panel- 
Projected positions in the plane-of-the-sky (x'-y' plane). Upper right and lower left panel: Declination (y') and Right Ascension (x') 



against line-of-sight velocity. Similarly to the observations, we apply a column density threshold of iVga 
distribution. 



10 cm to the simulated gas 



grid niodeling of star formation or whether the previous 
dynamical models (e.g. initial conditions) were not ac- 
curate enough. 

A first s imulation of the Antennae galaxies was pre- 
sented by iToomre fc Toomr^ ()1972D . reproducing the 
correct trend s in the morphology of the tidal tails. 
I Barnes! (|1988| ) repeated the analysis with a self-consistent 
multi-component mode l consisting o f a bu lge, disk and 
dark halo component. iMihos et al.l (|1993f ) included gas 
and star formation in their model and found the star 
formation to be concentrated in the nuclei of the disks, 
thus, not reproducing the overlap star formation. 

In this Letter, we present the first high-resolution 



merger simulation of NGC 4038/39 with cosmologically 
motivated progenitor disks galaxy models. We are able 
to match both the large-scale morphology and the line- 
of-sight kinematics, as well as important key aspects of 
the distribution and ages of newly-formed stars at the 
center of the Antennae, being a direct consequence of 
the improved merger orbit. 

2. SIMULATIONS 

The simulation presented here is the best-fitting model 
of a larger parameter stud y and was performed using 
Gadget 2 (|Springelll2005[ ). We include primordial ra- 
diative cooling and a local extra-galactic UV background. 
Star formation and associated SNII feedback are mod- 
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TABLE 1 

Model parameters of the best-fit merger configuration 
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"Mass in lO^" Mq 

'Disk and bulge lengths (r^iski ''bulge) a^nd disk scale height (zq) 
are given in kpc 
'^Halo concentration parameter 
''Halo spin parameter 
'^Maximum rotational velocity in kms~^ 

eled following the sub-grid multi-phase prescription of 
iSpringel fc Hernguistl ()2003l) . but excluding supernovae- 
driven galactic winds. For densities n > ncrit = 
0.128 cm~'^ the ISM is treated as a two-phase medium 
with cold clouds embedded in pressure equilibrium in a 
hot ambient medium. We deploy a fiducial set of param- 
eters governing the multi-phase feedback model result- 
ing in a star formation rate (SFR) of ~ 1 M0 yr~^ for a 
Milky Way-type galaxy. We adopt a softened equation 
of state (EQS) with gEQS = 0.5, where the parameter 
9EQS interpolates the star formation model between the 
full feedback model (<7eqs = 1-0) and an isothermal EQS 
with T = 10-* K ((jEQS = 0) (see iSoringel et all [20051 
for further details). The progenitor gala xies are set up 
in equ ilibrium according to the method of lSpringel et all 
(|2005| ) with a total virial mass of AUi, = 5.52 x 10" Mq 
for eac h galaxy. T he dark matter halos are constructed 
using a lHernguistl ([1990) density profile. They are popu- 
lated with exponential stellar disks comprising a constant 
disk mass fraction = 0.075 of the total virial mass and 
a stellar Hernquist bulge with a bulge mass fraction of 
mb — 0.025 (mb = l/3md). The gas mass fraction of the 
disk component is /g = 0.2 with the rest of the disk mass 
remaining in stars. T he disk and bulge scale len gths are 
determined using the iMo. Mao, k White! ()1998[ ) formal- 
ism. A summary of the most relevant model parameters 
is given in Table [TJ 

Each galaxy is realized with iVtot — 1-2 x 10^ par- 
ticles, i.e. 400,000 halo particles, 200,000 bulge parti- 
cles, 480,000 disk particles and 120,000 SPH particles. 
In order to avoid spurious two-body effects we ensured 
that all baryonic particles have the same mass and only 
one stellar particle is spawned per SPH particle. This 
yields a total baryon fraction of /bary — 10% with par- 
ticle masses for the baryonic components (bulge, disk, 
formed stars, and gas) of mbary — 6.9 x 10** M0 and 
wdm = 1-2 X 10^ M0 for the dark matter particles. 

The gravitational softening lengths are set to ebary = 
0.035 kpc for baryons, and edm = 0.15 kpc for 
dark matter particles, scaled according to eDM — 

ebary (™DM /wbary ) • 

We adopt an initially nearly-parabolic, prograde orbit 
geometry (the orbital plane lies in the x-y plane) with a 
pericenter distance of Tp — rd,4038 + ''d,4039 = 10.4 kpc 



and an initial separation of rs^p — J'vir = 168 kpc. For 
the orientation of the progenitor disks we found the best 
match to the Antennae system with inclinations «4038 = 
60°, 14039 = 60° an d arguments of pericente r W4038 = 
30°, W4039 = 60° (see lToomre fc Toomrel[l97l . 

3. RESULTS 
3.1. The morphological and kinematical match 

We determine the time when the simulation best 
matches the Antennae together with the viewing angles 
{9,'ip,4') which specify a series of subsequent rotations 
around the x-, y-, and z-axis. In the further analysis we 
will use the rotated 3D position-velocity subspace, i.e. 
the plane-of-the-sky (x'-y' plane) and the line-of-sight 
velocity vios, for comparison with the observations. Fi- 
nally, we apply a dis tance scale (£) relativ e to a fiducial 
distance of 22 Mpc (|Schweizer et al.l[2008[ )^ and assume 
a systemic helio-centric velocity of 1630 km s~^ to fit the 
observational data to the physical scales in the simula- 
tion. 

We find our best match to the observed large- and 
small-scale properties of the system with viewing an- 
gles of (93,69,253.5) and jC — 1.4, yielding a distance 
oi D ^ 30.8 Mpc to the system. The "best fit" {t = 
1.24 Gyr after beginning of the simulation) is reached 
only ~ 40 Myr after the second encounter {t = 1.20 Gyr), 
and approximately 50 Myr before the final merging of the 
galaxy centers (t — 1.29 Gyr). From our larger param- 
eter study we found this exact timing to be a manda- 
tory requirement for reproducing the overlap starburst. 
The first close passage of the two progenitor disk galax- 
ies occurred ^ 600 Myr ago which is considerably longer 
ago than ^ 200 — 400 Myr as sugg ested in earlier mod- 
els ()Barneslll98"l: iMihos et al']|1993l) and in much better 
agreement with observed 'intermediate-age' star clusters 
(^ 500 — 600 Myr). In FiglT] we show three large-scale 
projections of the PV cube of our simulated gas particles 
(NGC 4038: blue, NGC 4039: red) togeth er with a direct 
compa rison to Hi observations (yellow) bv lHibbard et al.) 
()2001|) . The Hi gas phase is used here as a tracer for 
the smooth underlying morphological and kinematical 
structure of the gas in the Antennae and we apply, sim- 
ilarly to the Hi observations a column density threshold 
of < 1 X 10^°cm~^ in the simulation. The top 

left panel displays the plane-of-the-sky projection, while 
in the top right and bottom left panels we show two 
orthogonal position-velocity profiles. Declination versus 
'^los (upper right) and vios versus Right Ascension (lower 
left). The simulation matches the morphology and kine- 
matics of the observed system very closely, especially for 
the southern arm, including the prominent kink in the 
velocities at the tip of the tidal arm (see FigdJ upper 
right and lower left panels). Due to the different initial 
orientations of the progenitor disks, the gas distribution 
in the northern arm is more diffuse than in the south- 
ern arm. The assumed column density cut-off therefore 

^ Note that the distance to the Antennae is a matter of recent de- 
bate. The systemic recession velocity yields a distance of 19.2 Mpc 
(assuming Hq = 75 kms~^ Mpc ~ ) while photometry of the tip 
of the red giant branch suggests a much shorter distance of only 
13.3 Mpc (Saviane et al. 2008^. Recently, Schwoizcr et al. (20081) 
have used three independent methods to determine a distance of 
22 ± 3 Mpc . 
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Fig. 2. — Line-of-sight velocity fields inside 18 kpc of the simulated and observed central disks. Isovelocity contours are drawn at 10 km s~^ 
intervals ranging from — ISOkms^^ to 150kms~^. Left: density-weighted velocity map. Right: intensity-weighted Hi velocity field of the 
high-resolution data cube. A column density threshold is applied as in Fig[T] 



results in a similar characteristic stubby geometry as ob- 
served (FiglU upper left). 

A closeup of the simulated and observed line-of-sight 
gas velocity fields in the central 18 kpc of NGC 4038/39 
is shown in the left and right panels of FiglJl Gas parti- 
cles are binned on a SPH-kernel weighted 256'^ grid and 
summed up along the line-of-si ght to prodiice a d ensity- 
weighted velocity map (see Hi bbard et al.l l2001t ) . The 
grid is smoothed with the observed beam profile and dis- 
played using the same p rojected pixel si z es (A ra = 2.64" 
and AdccI = 2.5") as in lHibbard et al.l (|2001[ ). We over- 
lay isovelocity contours spaced by lOkms"^ and apply 
the same column density threshold as in FiglT] The sim- 
ulation agrees well with the observed velocity field of the 
disk of NGC 4038. The northern part is approaching and 
the southern part is receding at similar velocities. Simi- 
larly, the simulated disk of NGC 4039 is approaching in 
the northern part and receding in the southern part like 
in the observations. In the simulation we have signifi- 
cantly more gas in the overlap region and the southern 
disk than in the observed Hi velocity field. This is due 
to the fact that we do not distinguish between molecular 
and atomic gas in our simulation whereas most of the 
gas in the centr al regions of the A ntennae seems to be in 
molecular form ()Gao et al.ll2001l ). 

3.2. The recent starburst 

In FiglU we show a color-coded map of the total gas 
surface density in the central 18 kpc of the simulation 
(upper panel). The nuclei of the progenitor disks are still 
distinct and connected by a bridge of high density gas. In 
the lower panel of Figl3] we show the corresponding iso- 
density contours and overplot in color all stellar particles 
(with individual masses of Wstar = 6.9 x lO^M©) formed 
in the last t < 15Myr (blue), 15Myr < t < 50Myr 
(green), and 50Myr < t < 100 Myr (red). In regions of 
currently high gas densities the very young stars (blue) 



form predominantly at the centers, in the overlap region, 
as well as in the spiral features around the disks si milar to 
the o bserved system (jWhitmore et al.i 1999 : Wang et al.l 
[200l . save the fact that the star formation in the centers 
seems to be much more pronounced in our simulation (see 
below). However, the overlap region is almost devoid of 
stars older than 50 Myr (red) indicating that the overlap 
starburst is a very recent phenomenon. Simulating the 
system further in time we find that the total duration of 
t he off-center s tarburst is no longer than « 20 Myr. 

iBrandl et al"] (2009) derived SFRs in the nuclei of NGC 
4038 (O.63M0yr-i) and NGC 4039 (0.33 Mq yr"!), and 
a total SFR of 5.4MQyr~-'^ for 5 infrared peaks in the 
overlap region. Comparing these values to simulated 
SFRs of 2.9 Mq yr-i (NGC 4038) and 2.8 Mq yr"! (NGC 
4039) in the galactic nuclei (defined as the central kpc), 
together with l.OMQyr"^ in the overlap region, we find 
that our simulation still falls short of producing the most 
intense starburst in the overlap compared to only mod- 
est star formation in the nuclei. We note, however, that 
we find a ratio (SFRovcriap/SFRnucici) of a factor of ^ 60 
(!) higher than reported in a previous Antennae model 
(Mihos et al. 1993). The total SFR of 8.1M0yr-i mea- 
sured from the SPH particles is in good agreement with 
the range of observ ed values between 5 — 20 Mq yr~^ (e.g. 
IZhang et al']l2001l ). 

In FigU] we plot the formation rate of stellar parti- 
cles within 18 kpc against their age (solid line). We 
find a significant increase of the SFR after the first and 
in particular after the second pericenter (dotted and 
dashed horizontal lines). Assuming the simulated SFR 
to be directly proportional to the cluster formation rate 
we compare the simulated SFR to obs ervations of the 
age distribution of yo ung star clusters ()Fall et al.l 120051 : 
IWhitmore et al.ll2007| ). using the same time binning of 
0.5 dex in log{T[yr\). We find that our simulated data 
are in very good agreement exhibiting a similarly good 
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match to the observe d cluster formation rat e as found by 
iBastian et alJ (|2009[ ) who compared to the iMihos et alJ 
(|1993() Antennae model. This model predicted a nearly 
constant formatio n rate for ages r < 100 Myr. However, 
in contrast to the ' Mihos et al.l (jl993) model, we find an 
additional significant increase in the formation rate of 
young stellar populations at ages r < 10 Myr induced by 
the recent second encounter. Despite the increase, the 
predicted formation rates of young clusters are still an 
order of magnitude lower than observed. Further inves- 
tigations have to show whether this discrepancy origi- 
nates from still uncertain details of the star formation 
model or explicit effects of the early disruption and evo- 
lution of massive clusters ("infant mortal ity", see e.g. 
IWhitmore et al.l l2007t IBastian et al] |2009( ) which were 
not included in our model. 

In this simulation, and other simulations in our param- 
eter study with similar central properties, we only find 
prominent star formation in the overlap region for a very 
short period of time after the second encounter, lasting 
for only ^ 20 Myr. In addition, stellar feedback is re- 
quired to prevent the rapid consumption of gas by star 
formation at earlier times. In a comparison run with- 
out stellar feedback most of the gas is consumed effi- 
ciently after the first pericenter and not enough gas is 
left over to form the overlap starburst after the second 
encounter. Thus, a central conclusion of our study is 
that the strong localized, off-center starbursts observed 
in the overlap region stems from a short-lived transient 
phase in the merging process associated with the recent 
second encounter. 

4. DISCUSSION 

The new numerical model for the Antennae galaxies 
presented in this Letter improves on previous models in 
several key aspects. We find an excellent morphological 
and kinematical match to the observed large-scale mor- 
phology and Hi velocity fields (|Hibbard et al.ll200ll ). In 
addition, our model produces a fair morphological and 
kinematical representation of the observed central re- 
gion. A strong off-center starburst naturally develops 
in the simulation - in good qualitative and quantitative 
agreement with the observed extra- nuclear star-f orming 
sites (e.g. [Mirabel et al.lll998t IWang et al.ll2004[ ). This 
is a direct consequence of our improved merger orbit. 
All previous studies using traditional orbits failed to re - 
produce the overlap starburst (see e.g. lKarl et al.|[2"008l ). 
The exact timing after the second encounter shortly be- 
fore the final merger ensures that the galaxies are close 
enough for the efficient tidally-induced formation of the 
overlap region. The formation of the extra-nuclear star- 
burst is likely to be supported by compressive tidal forces 
which can dominate the overlap region in Antennae-like 
galaxy mer gers during close encounters (jRenaud et al.l 
120081 120091 ). Energetic feedback from supernovae pre- 
vents the depletion of gas by star formation at earlier 
merger stages and ensures that by the time of the second 
encounter enough gas is left over to fuel the starburst. 
Simulating the system with an identical orbit, but now 
employing an isothermal EQS without feedback from su- 
pernovae (gEQS = 0) resulted in most of the gas being de- 
pleted by star formation at earlier phases of the merger, 
i.e. during the first encounter. 

Our model predicts that the observed off-center star- 
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Fig. 3. — Top: Gas surface density in the central 18 kpc of the 
simulation. There are clear concentrations of gas in at the two 
centers of the galaxies and the overlap region (red contours). Bot- 
tom: Recently formed stellar particles color-coded by their ages. 
The youngest stars (blue: r < 15 Myr) have formed predominantly 
in the overlap region and the centers associated with the peaks in 
the gas surface density, and tidal features around the disks (see 
upper panel). Older stars (green: 15 Myr < t < 50 Myr; red: 
50 Myr < r < 100 Myr) have formed throughout the galactic disks 
and the tidal arcs. 

burst is a transient feature with a very short lifetime (« 
20 Myr) compared to the full merger process (~ 650 Myr 
from first encounter to final merger). This fact serves as 
a plausible explanation for why such features ar e rarely 
observed in interacting galaxies (jXu et al.ll2000f ). How- 
ever, the observed puzzling gas concentration between 
the two nuclei of th e NGC 6240 merger system might be 
of a similar origin (iTacconi et al.l 119991 Engcl m prep. 
2010) suggesting that the Antennae overlap region, al- 
though rare, is not a unique feature. 
In addition, our improved model can serve as a solid 
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Fig. 4. — Formation rate of stellar particles versus age for our 
simulation (stars and solid line). Vertical lines indicate the time of 
first (dotted) and se cond (dashed) pericente r. The observed cluster 
formation rate from lWhitmore et al.l 120071) is given as filled circles. 

basis and testbed for further theoretical studies of the 
enigmatic interacting NGC 4038/39 system. For exam- 
ple, the overlap region in the Antennae is dominated by 
molecular gas, which we do not model in the simulation 
presented here. Given that we now have a dynamically 
viable method for forming the overlap region, detailed in- 
vestigations of the molecular gas formation process can 
be undertaken using improved theoretical model s (e.g. 
[Robertson fc Kravtsov..2008; .Pclupcssv fc PapadopoulosI 



l2009f ). In a first application using this new orbital con- 
figuration we have been able to qualitatively and quanti- 
tatively reproduce the magnetic field m orphology of the 
Antennae galaxies ([Kotarba et al.|[2009[ ). 

Finally, accurate modeling of nearby interacting sys- 
tems also provide unique insights into the merger dy- 
namics and timing of observed merger systems. The An- 
tennae galaxies are traditionally in the first place in the 
classical Toomre sequence which or ders galaxies a ccord- 
ing to their apparent merger stage ()Toomrelll977l) with 
the Mice (NGC 4676) being b etween their first and sec- 
ond pericenter (jBarnesI l2004[ ) and thus in second place 
behind the Antennae. According to our proposed model 
the Antennae galaxies are in a later merger phase, after 
the second pericenter. As a consequence the Antennae 
would lose their first place, and thus requiring a revision 
of the classical Toomre sequence. 
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